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We consider analysis targets at the International Linear Collider in which only a single 
photon can be observed. For such processes, we have developed a method which uses 
likelihood distributions using the full event information (photon energy and angle). 
The method was applied to a search for neutralino pair production with a photon from 
initial state radiation (ISR) in the case of supergravity in which the neutralino is the 
lightest supersymmetric particle. We determine the cross section required to observe 
the neutralino pair production with ISR as a function of the neutralino mass in the 
range of 100 to 250 GeV. 

1 Introduction 

Weakly interacting massive particles (WIMPs) are important study targets in current and 
future collider experiments since the WIMPs provide a solution to the dark matter problem. 
In this study, we investigate the prospects of discovering the direct pair production of WIMPs 
at the International Linear Collider (ILC) using the energy and angle spectra of the photons 
from initial state radiation (ISR) P]. We make the assumption that, among the particles 
beyond the Standard Model, the WIMPs are the only kinematically accessible at the centcr- 
of-mass (CM) energy of 500 GeV proposed for the ILC. Similar past searches using single 
photons have been performed by e.g. [2l|3]. 

Prospects of such measurements have been discussed in many literature [H [Sj |6l [7j and, 
in particular, at the ILC [8| which uses a cut-and-count approach. Our method uses the 
two-dimensional energy and angle distributions of the ISR photon and thus exploits the full 
information of the event. We apply this method to an example case of minimal supergravity 
(mSUGRA) [9l[Tni[ll] in which the neutralino is the lightest supersymmetric particle (LSP). 
The mSUGRA parameters have been chosen so that the LSP becomes bino-like; under this 
assumption, the dominant contribution to the neutralino pair production becomes the t- 
channel diagram shown in Figure [T] We perform a likelihood analysis based on the photon 
distributions and derive cross-section limits. We assume an integrated luminosity of 500 fb~^ 
at = 500 GeV. 

2 Signal and background samples 

Signal events have been generated using WHIZARD [H]. A GEANT4-based [13] simulation 
of the detector response was carried out by Mokka [M] using the International Large Detector 
(ILD) Concept [15]. For the background studies, we used the full Standard Model samples 
which were generated for the ILD Letter of Intent [15] . 

Signal samples were generated with two different LSP masses by varying the mSUGRA 
parameter M1/2, resulting in LSP masses of approximately 150 and 200 GeV. The other 
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Figure 1: Tree-level diagrams for neutralino pair production with initial state radiation: 
s-channel (left) and t-channel (right). 

mSUGRA parameters have been fixed to the following values: itiq = 300 GcV, = 0, 
tan/3 = 5 and sgn(/Li) = +1. The mass spectra have been computed using SOFTSUSY [TC] : 
the relevant values are summarized in Table [T] The default cuts of WHIZARD have been 





Sample A 


Sample B 


m^o (LSP) 


151 GeV 


205 GcV 




395 GeV 


453 GeV 




333 GeV 


354 GeV 


Ml/2 


375 GeV 


500 GcV 



Table 1: Mass spectra of signal samples. 



modified to allow the generation of photons with energy greater than 0.1 GeV and | cos 6*1 < 
0.995 where the angle 9 is taken to be between the direction of the photon momentum and the 
beam axis. The signal samples have been generated separately for different polarizations for 
the electron and positron beams. For the case of m^o = 150 GeV, the cross sections restricted 
to the energy and angle ranges used for the event generation are a{e~j^e\^ — ^ 7XiXi) = 56 fb 
and a{ej^e^ — > 7X1X1) = 1.9 fb; for the case of 200 GeV, the respective cross sections 
are 18 fb and 0.55 fb. The difference in the cross section between the two polarization 
combinations arises primarily due to the weak hyperchargc of the chiral fields and, to a 
lesser extent, the difference in the scalar electron masses. The two samples were mixed with 
appropriate weights to yield a sample with the beam polarizations {Pe~ , Pe+) = (+0.8, —0.3). 
We list the dominant background processes after the event selection in Table [2J 



Process 


Cross section 


77 t+r 

e^e~ — >■ ^vivi 


1.1 X 10^ fb 
7.8 X 10^ fb 



Table 2: Dominant background processes: the cross sections in the table are calculated 
after the event selection as described in Section [3l The lepton generations (£ = e,/x,r) are 
summed in the cross section numbers. 
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3 Event reconstruction 



We rely on the Pandora Particle Flow Algorithm [T^ for the reconstruction of calorimeter 
clusters and charged particles. A single highly energetic photon often results into multiple 
clusters. In order to recover the split clusters, we merge the clusters which lie within a cone 
of an opening angle of 1.5 degree around the most energetic cluster. For the event selection, 
we require that there are no charged particles in the event. In addition, we require that 
there are no other clusters aside from the merged cluster. The merged cluster is required 
to have an energy greater than 0.5 GeV. For the particle identification, we take the energy 
deposits in the electromagnetic calorimeter i?ECAL and the hadronic calorimeter E-ucAC to 
compute the energy fraction £^ECAL/(i^ECAL + -Ehcal) which is required to be greater than 
0.9. The reconstruction efficiency of signal events as a function of energy and angle is shown 
in Figure [2] The overall signal efficiency is 94% after the event selection. 
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Figure 2: Reconstruction efficiency as a function of angle and energy after the photon cluster 
merging. 



4 Likelihood analysis 

We perform a binned likelihood analysis using the two dimensional distribution of the photon 
energy (logarithmic scale) and angle. We choose a binning such that the number of bins is 
N = 96, with the binning widths shown in Figure [31 The choice of log E instead of E is 
motivated by the fact that the ISR photon distribution diverges at low energy. The angle is 
folded around 9 = 7r/2, to exploit the symmetry of the distribution. The bin size is adjusted 
to ensure the presence of sufficient signal events in each bin. The log likelihood ratio is 
defined as follows: 

ln£(ni,n2,...,nAr) = 2^1n — — (1) 

where n,; is the observed number of events in the i-th bin, bi is the expected number of 
background events in the i-th bin, Si is the expected number of signal events in the i-th bin, 
and P(n; A) = A"e-^/n! is the Poisson probability density function with the expectation 
value of A. The m-^o = 150 GeV sample was used to build the template signal shape; the 
same signal shape was used to build likelihood functions even for different masses in the 
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signal <m=1SGGeV) 2(1 distribution ot the plioton Signal (m-20a6aV) 2d distribution of tha photon Background 2d distribution nf the phaton 




log^^ (Energy/GeV) log^^ (Energy/GeV) log^^ (Energy/GeV) 

Figure 3: Two-dimensional distributions for signal events with m^o = 150 GeV (left), 
m^o = 200 GeV (center), and background events (right). The horizontal axis is the photon 
energy in logarithmic scale. The vertical axis is the photon angle in radians, folded around 
7r/2. The binning is shown as red lines; the lower-left corner is cut off to make the binning 
consistent with the cuts in the background sample. 

subsequent analysis, since the true mass is not known a priori. The likelihood distributions 
are obtained by performing Monte-Carlo test experiments whereby the observed number of 
events in each bin is fluctuated using Poissonian random numbers. The resulting likelihood 
distributions for background plus signal events (Ai = bi + Si) are compared with that for 
background events only (A^ = bi) to compute the probability for observing the WIMP 
signal, as shown in Figure |4l From this distribution, we calculate the expected significance 
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Figure 4: Distribution of log-likelihood values for background only (red), background and 
signal with m^o = 200 GeV (green), and background and signal with m^o =150 GeV (blue). 
The vertical axis is the probability density for a given log-likelihood value. 

for each neutralino mass. We estimate the significance by computing the probability that 
the backgrounds fluctuate to deviate more than the median of the signal likelihood. In the 
case of m-^o = 200 GeV, the estimated significance is S.lcr; in the case of m^o = 150 GeV, 
the estimated significance is well above 5a. The difference in the significance between the 
two cases is mostly due to the difference in the cross section. 
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In the case of a more general WIMP scenario, the dependence of the signal cross section 
on the WIMP mass is expected to be considerably different. In order to take this effect 
into account, we performed a scan of likelihood values by varying the WIMP mass and the 
cross section independently. The best signal separation power can be attained when the 
correct signal template function is used; assuming this is the case, we keep our mSUGRA 
signal template. Additional signal samples were generated using 15 different WIMP mass 
points ranging from 100 GeV to 240 GeV by varying Mi/2- For this study, we use the beam 
polarizations of (Pg- , Pe+) = (+1-0, —1.0). We obtain the cross section limit, defined as the 
cross section which gives a 50% probability of 5ct observation, as a function of the WIMP 
mass. 
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Figure 5: Cross section limits for WIMP observation. The vertical axis is the cross-section 
limit of the polarization combination (Pg-,Pe+) = (+1.0,-1.0) integrated over the kinc- 
matically allowed region in this analysis. 

For a given mass, we adjusted the cross section until it yielded the desired probability 
of observing the signal. We performed 10^ Monte-Carlo test experiments at each step of 
testing the cross section. A binary search was used to find the cross section limit efficiently. 
The cross section limit as a function of the WIMP mass is plotted in Figure [SJ If the 
cross section is at least 14 fb, the WIMP can be observed. Around the mass of 200 GeV, 
the required cross section increases since the likelihood function is constructed using the 
template signal shape for m^o = 150 GeV. As the WIMP mass approaches 250 GeV, the 
required cross section becomes smaller since the photon energy distribution peaks at the low 
energy region, thereby allowing it to be discriminated against the background shape. Note, 
however, that if one fixes the model and let the mass approach 250 GeV, the cross section 
decreases faster than the cross-section limit. Thus, this result is consistent with the fact 
that the heavier neutralino is more difficult to observe. 

5 Conclusions 

We have looked at the mSUGRA neutralino search as an example of a WIMP search with 
single photons using the two-dimensional distribution of the ISR photon energy and the 
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angle. We compared the likelihood distributions and calculated the probability to achieve 
5a observation. In our mSUGRA model, the m^o = 150 GeV case can be easily observed, 
while the m^^o = 200 GeV case has the significance of S.lcr. We obtained cross-section limits 
as a function of the neutralino mass. WIMPs created through the ^-channel with the cross 
section of over 15 fb can be observed with the 5cr significance at the ILC with y/s = 500 
GeV. In principle, this method can be applied to a broader class of WIMP searches with 
ISR photons. 
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